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The thermodynamics of low dimensional spin-1/2 Heisenberg ferromagnets (HFM) in an external 
magnetic field is investigated within a second-order two-time Green function formalism in the wide 
temperature and field range. A crucial point of the proposed scheme is a proper account of the 
analytical properties for the approximate transverse commutator Green function obtained as a 
result of the decoupling procedure. A good quantitative description of the correlation functions, 
magnetization, susceptibility, and heat capacity of the HFM on a chain, square and triangular lattices 
is found for both infinite and finite-sized systems. The dependences of the thermodynamic functions 
of 2D HFM on the cluster size are studied. The obtained results agree well with the corresponding 
data found by Bethe ansatz, exact diagonalization, high temperature series expansions, and quantum 
Monte Carlo simulations. 

PACS numbers: 75.40.Cx 



I. INTRODUCTION 

Quantum Heisenberg model with ferromagnetic ex- 
change is used extensively to interpret the ther- 
modynamic and magnetic properties of low dimen- 
sional (ID and 2D) physical systems. Examples 
of quasi-one-dimcnsional ferromagnets whose proper- 
ties can be explained within t he Heis enberg model 
are organic p-NPNN compound ;] 1 ! 2 ! 3 ! 4 ^ and cuprates 
TMCuCP 3 The ferroma gnctic insulators such as 
K 2 CuF 4 , Cs 2 CuF 4 , La 2 B aCuQ 5 R^CrCl^^ and 
quantum Hall ferromagnet d 11 * 12 -^ provide examples of 
the Heisenberg system on a square lattice. A unique 
example of a spin-1/2 magnet on a triangular lattice is 
3 He bilayers adsorbed on graphite.E^ 1B ' 17 ' 18 ' 19 ' 2U ' 21 J At 
high coverages the second layer proved to be a solid fer- 
romagnet whose thermodynamics ca n be describ ed by 
HFM with a high degree of accuracy ft 1 | 22 | 23 | 24 | ^l Nowa- 
days, considerable study is being given to 3 He monolayers 
on 4 He-preplated graphite substrates. In these systems 
under high enough pressure a solid 3 Hc monolayer with 
ferromagnetic exchange is formedP^ZI 

Experimental research of the aforementioned low di- 
mensional magnets is carried out intensively, in partic- 
ular, in the presence of an external magnetic field. To 
interpret the experimental data it is necessary to develop 
a quantitative description of the HFM thermodynamics 
at arbitrary magnetic fields and temperatures. The two- 
time Green function formalism is quite appropriate for 
this purpose. The method based on one or another de- 
coupling scheme for higher Green functions results in a 
closed se t of self-co nsistent equations for thermodynamic 
averages! 28 * 29 * 30 * 3 -^ Random phase approximation (RPA) 
is the simplest variant of such scheme with decoupling 
at the first steppe Being applied to the low dimensional 
systems it gives satisfactory results at high fields whereas 
at low and intermediate fields RPA describes the thermo- 
dynamics only on a qualitative level. 

A quantitative description of ID and 2D spin systems 



can be obtained within a more complicated scheme orig- 
inally proposed in Ref. [33] for ID HFM in zero mag- 
netic field. The scheme is based on the decoupling of 
higher Green functions at the second step with intro- 
ducing the vertex parameters to be found. A proper 
choice of the vertex parameters makes it possible to re- 
tain some relations that must hold true at the exact solv- 
ing of the problem. As a result, the theory is built in 
terms of the correlation functions and the vertex param- 
eters obeying the self-consistent set of equations. In Refs. 
1331341351361371381 a single vertex parameter was chosen so 
as to satisfy the sum rule. One vertex parameter turned 
out to be quite enough to describe quantitatively the 
thermodynamics of ID and 2D (on a square and trian- 
gular lattices) ferromagnets in zero field. 

When employing the above-mentioned scheme to the 
spin systems in an external magnetic field, along with 
the correlators we have to determine at least three ad- 
ditional functions of temperature and field: two vertex 
parameters and magnetization. To do this we need three 
relations two of which arc quite evident from the proper- 
ties of the spin-1/2 operator 

((S z ) 2 ) = \, {S Z ) = \-{S-S+), (1) 

where angular brackets denote thermodynamic averag- 
ing. The choice of the third condition is not so apparent. 
In Ref. [35] a second-order Green function scheme was ap- 
plied to HFM chain and HFM on a square lattice. As the 
third condition the authors of Ref. [39]used the exact rep- 
resentation of t he int ernal energy through the transverse 
Green functionPHUl 

The aim of the present work is to calculate the 
thermodynamic functions of HFM on a triangular lat- 
tice in an external magnetic field using a second-order 
Green functio n formalism. As compared to a chain and 
square lattice] 12 ! 39 " * 41 ! 42 ! HFM on a triangular lattice 
in a magnetic field is much less investigated. For this 
case high temperature series expansion (HTSE)J^ low- 
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temperature asymptotics for the magne tizati on calcu- 
lated within the spin wave approximation E4J45I tempera- 
ture dependences of the magnetization found by quantum 
Monte Carlo simulations (QMC) on a 16 x 16 cluster,^ 
and some results obtained by the renormalization group 
technique^ (RGT) are known. However, none of these 
approaches gives a complete description of the thermo- 
dynamics in the whole temperature and field ranges, and 
a second-order Green function method is expected to fill 
the gap in our knowledge. 

In the present work we show that this method is more 
effective when the conditions determining the magneti- 
zation and vertex parameters result from the fundamen- 
tal principles. Clearly the relations ([!]) are just of this 
kind. It is equally important to retain the analytical 



properties 
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of the Green functions in the approx- 



imate approach. Such a requirement for the transverse 
commutator Green function provides a basis for the third 
condition in our theory. Note also that by appropriate 
choice of variables the set of self-consistent equations for 
the correlators, vertex parameters, and magnetization 
can be written in a universal form suited not only for 
a triangular lattice but also for a chain and square lat- 
tice. A good agreement of our results obtained for the 
three types of lattices with the corresponding data avail- 
able from literature confirms the efficiency of the used 
scheme. 

In Sec. II, the statement of the problem is formulated 
and the self-consistent set of equations for the correla- 
tors, magnetization, and vertex parameters is derived. 
In Sec. Ill, the proposed scheme is applied to HFM on 
a chain, square and triangular lattices. The obtained re- 
sults are compared with the corresponding data found 
within other methods. Some concluding remarks are 
made in Sec. IV. 



II. STATEMENT OF THE PROBLEM 



The Hamiltonian of the system is given by 



J 



(2) 



where a± and a z are the vertex parameters. 

After a number of manipulations we finally obtain for 
the time-space Fourier component ((S'^|S'i k )) aJ 



Jci7o i-r k 



4tt w 2 -(w£) 2 ' 



(4) 



where 



K) 2 - - r k ) [A z + 70^(1 - r k )] , (5) 



A z = 1 + c 2 - (70 + l)ci. 



(G) 



Here the following correlation functions have been intro- 
duced 



c x = 2(SfS f ° s ), c 2 = 2^ '(S? +s S f ° s ,), ci,2 = a z ci, 2 . 

The primed sum indicates that the term with S = S is 
omitted in it. The structure factor r k is defined as 



r k = — y^exp(zk^), 
To * 



(8) 



where the coordination number 70 is equal to 2 for a 
chain, 4 for a square lattice, and 6 for a triangular lattice. 
Fourier transform ((S^.\SZ^))u can be written as 



Here 



1=1,2 '. k 



(9) 



(10) 



(^k) 2 = ^?(i-r k ) a -,,/Mi-r k i . in, 



t,s 



l + b 2 -(70 + 



(12) 



where Sf is the spin-half operator at site f , d is a vector 
connecting nearest neighbors, J > is an exchange inte- 
gral, h = 2{iB, fi is the magnetic moment of a particle, 
B is an external magnetic field. 

To calculate spin-spin correlators, it is necessary to 
find two retarded commutator single-particle Green func- 
tions: ({SI\SI,}), ((S?\S f ,' T )) (a = ±). We write down 
equations of motion for these two functions and make the 
decoupling of the higher Green functions on the second 
step according to the scheme proposed in Ref. [39] 



S^SjSf — a±(Sf Sj)S[ , S^S, a S* — a z (SfS, a )S[ , 



(3) 



4 <T — 



a \ S H ^~~L U^k), 



(13) 



where (S z ) is the magnetization. Due to the presence of 
the external magnetic field, (S z ) is nonzero at any finite 
temperature. 



The correlation functions entering Eqs. (Ill— (131 are 
defined by 

b t = ^^, k = a x b l , 1 = 1,2, (14) 



01 = A{SISI +S ), a 2 =4^ '(S' f+s St +5 ,). (15) 



The Green functions look formally the same for three 
above-mentioned types of lattices. Such universal form 
has been possible to obtain, because instead of the usual 
correlators describing correlations between spins which 
are two steps along the translation vector 5 apart, we 
use linear combinations ci and a 2 defined by Eq. Q and 
Eq. (15). The physical meaning of these combinations 
depends on the lattice type. For a chain, C2 and a 2 are 
the next nearest neighbor correlators. For a square lat- 
tice these combinations contain the correlation functions 
between the spin at site f and spins from the second and 
third coordination spheres, and for a triangular lattice 
these combinations in addition to the higher order cor- 
relators include also c\ and a\. For a chain and square 
lattice Green functions Q and (|9| coincide with those 



found in Ref. M 

Using the spectral relations^ we have 



z\2 



° 2 = ^ E i^Tl - 1) , 9k + 4 ( 7o - 1) (S z 



N 



Cl 



^E r ^ 



a±N 



where 



, k = ^(i-r k) coth(M 

u>i V 2 



^ = 2a ± (S z ) sinh (0h) - Jb l7 o (1 - T k ) sinh (0lj£) /ujj 
cosh (f3h) — cosh ((3uj^) 



N is the total number of sites, (3 = 1/T. Eq. (16 1 rep- 



resents the set of equations for the correlation functions 
q and ai. Along with these correlators, the set (16 1 con- 
tains the parameters a z , a±, and magnetization (S z ) to 
be also determined. 

The vertex parameters are chosen so as to satisfy the 
sum rules 

4((S|) 2 ) = 1, 2<SfS f - ff > = i + 2a(S 2 ), 
which using (|17[) can be written as 



Jci 



£<7 k + 4<S*> 2 = l, ^ = ^E 



Let us analyze Eq. (19 1. The frequency uo^_ has a 
maximum to^ax a t the edge of the Brillouin zone 



KuJ = 7o J 2 (A ± + 2 7o 6 



(21) 



Since the parameters Ax and b\ in Eq. (21 1 are functions 
of temperature, the frequency u>max depends on tempera- 
ture as well. It can be shown that w^ aa . decreases mono- 
tonically from u>^ ax = J 7 o at T — to oj^ lax — J^/Jo at 
00. At h/J < ^/tq Eq. (19 1 has a real solution for 



T 

any temperature. Substituting it into Eq. (20 1 we obtain 



,Pk- (18) the following expression for the magnetization 



Finally, in order to close the system (16 1, (18 1 we need 
one more equation. It c an be found from the following 
consideration. It is knowrP^Sl that a commutator Green 
function must not have any pole at uj — 0. Clearly Green 
function Q does not have such pole. A different situ- 
ation arises with Green function ([9]). When w = its 
denominator is equal to zero at k = ko with wave vector 
ko satisfying the equation 

h = u£ - (19) 

Thus, the numerator of Green function ^ must also van- 
ish at k = k , for otherwise this function would have a 
pole. From this condition we get the equation for (S z ) 



J&170 



(20) 



Note, that in calculating the anticommutator transverse 
Green function, the condition (20 1 appears automatically 



without any special assumptions (see also Ref. I5T) 



(S z ) = 




-A 



(22) 



In the field range ^/7o < h/ J < 70 the real solution of Eq. 
( 19 1 exists only at T < Tq (where Tq obe ys t he equation 



.(T ) = h). Finally, if h/J > 7o Eq. (19 1 has no real 



solutions at any temperature. 

It is natural to suppose that the expression ( 22 1 for 
the magnetization is valid at arbitrary h and T. This 
assumption provides continuity of (S z ) as a function of 
field and temperature. Eq. (p2| gives correct values of 



the magnetization at low and high fields for arbitrary 
temperatures and at T = for arbitrary fields. How- 
ever, the most important thing is that Eq. (22 1 provides 



correct analytical properties of the commutator Green 
function (|9| obtained within the approximate scheme. 



As a result, Eqs. ( 16 1, ( 18 1, and ( 22 1 represent a closed 



set of seven self-consistent equations for a\, a 2 , C\, C2, a z , 
ox, and (S z ). This set can be reduced to three equations 
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for hi, A 2 , and Aj_ 



2bi = a±ci 



Ja±ci 



Jci 



-E 

k 



k 

(i-r k ) ffk 



2\2 



2a 2 



;(i-r k )(i-7or k )9k. (23) 



The efficiency of the proposed scheme can be esti- 
mated, first, by comparing the obtained results with 
available from literature data found by alternative meth- 
ods and, second, with the help of inherent criteria exist- 
ing within the developed scheme itself. The first criterion 
implies that at T — ► oo the entropy (per site) of the sys- 
tem with spin 1/2 should tend to S(oo) = In 2 0.693. 

The second criterion follows from the relational con- 



necting the internal energy ( 25 I and the Green function 



that can be written as 



The values ci, a±, and (S z ) can be expressed through 
b\ 7 A z and Aj_ according to (16 1, (18) and (22). For a z 
with the help of Eqs. ^ and ^7 1 we have 



(70 + l)ci- c 2 



(24) 



It is easy to see that the replacement h — > —h changes 
the sign of the magnetization and does not change the 
correlation functions and vertex parameters. Owing to 
condition p2l, Eq. (Il7) and thereby Eqs. (fl6|, (pi) 



have no singularities, At h = the system (23 1 reduces 
to that found in Refs. 1331381 

The internal energy E per site is given by 



E 



J Jo 



(2 Cl + ai )-h(S z ). 



(25) 



^{2c 1+ai -l) + h({S z )-\ 



1 E 



, {e k + u,)Tm((S+\SZj u 
I duJ tf^Tl = °' (26) 



J 7o (l - T k )/2 + h. The relation (26) be 



where e, 

comes the identity with the exact Green function and 
correlators. This is not the case for the Green function 
and correlators found as a result of the decoupling proce- 
dure. Dividing Eq. (26 1 by J^q{S z )/A and substituting 
(J9j) in its left hand side, we get 



Jlo{S')N 



E 



J6170 (1 - T k ) sinh (ph) - 2(S z }lj£ sinh (pu£) 
cosh (ph) — cosh (/3w k ) 



R 



(27) 



The quantity R is a function of field and temperature. 
Due to Eq. (22) expression (27 1 is singularity-free. It is 



evident that the closer R to zero, the better the approxi- 
mation. Thus, the condition \R\ <C 1 can serve as another 
criterion of the approximation efficiency. Below, in dis- 
cussing the results we will calculate R(h,T) and check 
the fulfillment of this criterion for the proposed scheme. 



III. RESULTS AND DISCUSSION 



In the general case the set of equations ( 23 1 can be 
solved only numerically. In limiting cases analytical re- 
sults could be obtained. At T = the system (23 1 gives 
correct values for the sought quantities 



Note, that employing a similar approach for HFM on 
a chain and square lattice (below we refer to it as Green 
function approximate method (GFAM)) the authors of 
Ref. EH1 instead of Eq. ^ used the condition R = 
as one of the equations in the self-consistent set of equa- 
tions. In the following we will compare the thermody- 
namic functions calculated within our scheme with the 
results found in Ref. 1391 



(S z 



C12 = 0, a x 



a-2 = 7o 



(28) 



The same solution is also true for finite temperatures at 
h — » 00. 

In the high temperature limit (J, h <^ T) the system 



( 23 1 can be solved by expanding in 1 /T. Restricting our 
consideration to the second order in x — J /(AT) and y - 



5 



h/(2T) we obtain the following asymptotic expressions 



cT = x+ T ( 7o 2 - 670 + 4), 
cT = \ [(70 - 2) (70 - 4) - x(7 2 - 14 7 o + 20)] , 

r = c a 2 s + (70 - i)y 2 , 

T 

a_ = a 1 = 1 — — 



„as 1 2 



<S 2 r = f(l+7oz) 



(29) 



For the heat capacity C = dE/dT to the third order in 
x and y we get 



C* as = 



37oz 2 



1 + 2 ho - 6 7o + 4 ) 

y 2 (1 + 3 7o x) 



(30) 



The expressions for the magnetization and heat capac- 
ity coincide with those obtained by the direct high tem- 
perature series expansion.^ As it follows from Eq. ( 29 ) 
the field-dependent terms in expansions for q occur in 
the fourth or higher order in l/T. 

We will demonstrate the efficiency of the proposed 
scheme applying it to the ID HFM and HFM on a square 
lattice. The main attention will be paid to low fields, be- 
cause it is this region that is the most difficult for the 
adequate description within approximate methods. 




FIG. 1: Temperature dependences of the magnetization for 
ID HFM at h/J = 0.005 and 0.05 (from left to right). The in- 
finite system: present theory (solid), BA^ (■) and GFAIVP 1 
(dotted). The cluster: present theory (dashed) and EEp2l(o). 



A. One-dimensional Heisenberg model 

In this subsection we consider ID HFM. Figs. [T]- 
[3] demonstrate the magnetization, susceptibility \ — 
d(S z )/dh, and heat capacity vs temperature at low fields 
obtained within our approach. In calculating the ther- 
modynamic functions on clusters, we use the periodic 
boundary conditions. The corresponding dependences 
found in Ref. [39] by Bethe ansatz (BA) and GFAM for 
an infinite chain and by the exact diagonalization (ED) 
for a cluster of 16 sites are also shown for comparison. 

For the magnetization and susceptibility our method 
yields a good agreement with the exact results in the 
whole temperature range. The positions of the maxima 
in our curves for x(^) coincide perfectly with those found 
within the exact methods and only a small difference in 
the peak heights is observed. 

At low fields the exact methods indicate the depen- 
dence of the thermodynamic functions on the chain 
length. As it is seen from Figs, [I] [2j at h/J = 0.005 
the curves (S Z (T)) and x( T ) obtained by ED for the 
finite-sized chain differ substantially from those found by 
BA for the infinite system. Our method gives a proper 
description of this effect. At higher fields (Figs. [T| [3| 
where the ED and BA results coincide our dependences 
for N — 16 and N — > 00 also coincide and show a good 
fit to the exact data. 

As it can be seen from Fig. [3]there is a certain disagree- 
ment between the heat capacities obtained within the ex- 
act and approximate methods at low fields. This result 
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FIG. 2: Temperature dependences of the susceptibility for 
ID HFM at h/J = 0.005. The infinite system: present theory 
(solid), BA^B) and GFAM 33 (dotted). The cluster: present 
theory (dashed) and EEP 1 (o). 

is quite understandable. Indeed, a similar scheme^ ap- 
plied to ID HFM at h — gives in the low temperature 
region sufficiently different run of the heat capacity than 
the exact solution. Nevertheless, even at h/J — 0.1 our 
theory not only gives the correct position of the max- 
imum but also reproduces a specific bend in the curve 
C(T) at T/J - 0.3. The agreement between C(T) calcu- 



0.0 0.5 1.0 1.5 

T/J 



12 3 4 
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FIG. 3: Temperature dependences of the heat capacity for 
ID HFM at h/J = 0.1. The infinite system: present theory 
(solid), BA^B), GFAIVP (dotted). The cluster: EEP(o). 
The inset shows C(T) at low fields ft/J=0.001, 0.005, 0.01, 
0.025, 0.05, 0.075, 0.1 (from bottom to top). 



lated within our scheme and found by the exact methods 
becomes better with increase in field. The inset in Fig. 
3 illustrates a double-peak structure of the heat capacity 
that within our method is identified at < h/J < 0.045. 
A similar structure of C(T) at low fields was first ob- 
tained in Ref. |39] 

We calculated the entropy for ID HFM. It tuned out 
that the higher is the field the closer is the limiting value 
S(oo) to In 2. For example, at h/J = 0.05 the entropy 
is S(oo) ~ 0.631 and at h/J = 1 it is 5(oo) ~ 0.687. 
The quantity R(h, T) was also found. At low and high 
temperatures it is practically equal to zero, so that the 
Green function ^ and correlators calculated within our 
method may be considered as satisfying Eq. (26). 



At 



given field in the intermediate temperature range where 
the correlators vary rapidly the quantity R is at maxi- 
mum. At h/J — 0.05 the maximum value of R is ~ 0.027 
whereas at h/J — 1 it does not exceed 0.006. With in- 
crease in h/J the quantity R decreases and the difference 
between results obtained by the exact and approximate 
methods vanishes. 



B. HFM on a square lattice 



FIG. 4: Temperature dependences of the magnetization for 
HFM on a square lattice at h/J = 0.1 and 0.4 (from left to 
right). The infinite system: present theory (solid), QMCp^ 
(A) and GFA1VP 1 (dotted). The 4x4 cluster: present the- 
ory (dashed) and ECP ( q ). The inset shows (S z ) vs T/J 
at h/J=0.05, 0.1, 0.2, 0.32, 0.4 (solid) in comparison with 
QMCP (A) and HTSEP (dashed) (from left to right). 




Now we proceed to the HFM on a square lattice. Fig. 
[4] demonstrates the magnetization as a function of tem- 
perature at h/J =0.1 and 0.4 for 4x4 and 32 x 32 
square lattice clusters together with QMCP ED,^ and 
GFAM 39 data. The inset shows the magnetizations for 
the 32 x 32 cluster at h/J = 0.05, 0.1, 0.2, 0.32, 0.4 found 



FIG. 5: Temperature dependences of the susceptibility for 
HFM on a square lattice at h/J = 0.2, 0.4, 1.0 (from top to 
bottom). The infinite system: present theory (solid), HTSlP^ 
(dash-dotted), and GFAM^l (dotted). The 4 x 4 cluster: 
present theory (dashed) and EEp^ (o). 
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FIG. 6: Correlation functions a\ (a) and ci (b) for HFM on 
a square lattice at h/ J=0.1 and 0.4 (from left to rig ht). The 
infinite system: present theory (solid) and GFAM^* (dotted). 
The 4x4 cluster: present theory (dashed) and EEp3 (o). 



within the present method, QMC,E> and HTSEP The 
comparison between our results and HTSE is possible, 
because at these values of field the magnetizations for 
the 32 x 32 cluster proved^ to be identical to those for 
the infinite lattice. It is seen that the temperature de- 
pendences of (S z ) obtained within our approach are in 
good agreement with the exact results. The proposed 
scheme reproduces correctly the dependence of (S z ) on 
the size of the system as well. The difference between the 
present results and GFAM in Fig. |i](see also Figs. [l}{3]) 
testifies that Eq. (22 1 is more preferable as compared to 
the condition R = for a quantitative description of the 
low dimensional HFM at small fields. It is also evident, 



that the thermodynamic functions of the square lattice 
HFM are more sensitive to the choice of the condition 
for (S z ) than the thermodynamic functions for ID HFM. 
Naturally, such a choice is expected to be even more crit- 
ical for the lattices with larger coordination numbers (for 
example, a triangular lattice). 

Fig. [5] illustrates the susceptibility x(^) f° r the square 
lattice together with \{T) found by HTSEP EDP and 
GFAMPI 

The temperature dependences of the correlation func- 
tions ai and C\ at h/J = 0.1 and 0.4 calculated for the 
4x4 cluster as well as for the infinite lattice are presented 
in Fig. [6] The ED and GFAM results are added for com- 
parison. At low fields a clearly defined dependence on 
the size of the system is seen. Our results for the infinite 
lattice differ noticeably from GFAM. For the 4x4 clus- 
ter a good agreement with the dependences calculated by 
ED is observed. 

The limiting value of entropy for the HFM on the 
square lattice is S'(oo) = 0.651 at h/J = 0.05 and 
S(oo) — 0.684 at h/J = 1, which is very close to the 
exact value In 2. The maximum value of R is ~ 0.058 at 
h/J = 0.05 and ~ 0.014 at h/J = 1. 

Thus, the results of subsections A, B show, that the 
theory based on the correct accounting for the analyti- 
cal properties of Green functions gives an adequate de- 
scription of the thermodynamic functions for the systems 
under consideration in the wide field and temperature 
range. 



C. HFM on a triangular lattice 

In this subsection we consider HFM on a triangular 
lattice in an external magnetic field with peculiar at- 
tention concentrated on small and intermediate fields. 
Fig. [7] represents the temperature dependences of the 
magnetization at different values of h/J. It is seen that 
our results agree closely with HTSE 43 up to the point 
of HTSE applicability, whereas the RPA curves coincide 
with HTSE only at relatively high temperatures. In the 
intermediate temperature range the RPA results differ 
sufficiently from ours even at h/J— 1.5 reproducing the 
temperature behavior of (S z ) only qualitatively. We have 
also calculated the magnetization at low temperatures 
using the renormalization group technique. The results 
obtained by both our approaches are in good agreement. 

Fig. [8] illustrates the temperature behavior of the sus- 
ceptibility. Analysis shows that with increase in field the 
maximum in x(T) decreases and shifts to higher tem- 
peratures. At h/J > 0.1 the height of the maximum 
as a function of h/J with a great degree of accuracy is 
described by a power law 



'h 
" ' J 



0.1696, b = -0.8634. (31) 



Temperature dependences of the correlation functions 
a\ and c\ at different h/J are shown in Fig. [9] Beginning 
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FIG. 7: Temperature dependences of the magnetization for 
HFM on a triangular lattice at h/J = 0.1, 0.3, 0.5, 0.7, 1.0, 
1.5 (from left to right). The present theory (solid), HTSEP 1 
(dashed), RGT (thin lines), and RPA (dotted). 
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FIG. 9: Correlation functions a± (a) and ci (b) for HFM on 
a triangular lattice at h/J — 0.1, 0.3, 0.5, 0.7, 1.0, 1.5 (from 
left to right). The present theory (solid) and RPA (dotted). 
Thin lines correspond to h —0. 



FIG. 8: Temperature dependences of the susceptibility for 
HFM on a triangular lattice at h/J = 0.2, 0.4, 1.0 (from top 
to bottom). The present theory (solid), HTSEP 1 (dashed), 
and RPA (dotted). 



with T/J ~ 0.5 the RPA results differ from ours suffi- 
ciently. It is easy to verify that at T > J the correlator 
a\ calculated within RPA to the first approximation in 
J/T is negative and equal to — J/(4T). Thus, almost 
at all temperatures RPA fails to describe correctly the 
correlation functions and, hence, the energy and heat ca- 
pacity. 



the heat capacity C(T) in comparison with HTSE 43 and 
RPA. It is seen that our results are in good agreement 
with HTSE. With increase in field the position of the 
maximum in the curve C(T) shifts to higher tempera- 
tures and its value first increases rapidly and then 
decreases. A similar behavior occurs for the heat capac- 
ity maximum on a square lattice. Ma xim um values 
and vs field are illustrated in Fig. 11 Field de- 



Fig. 10 demonstrates the temperature dependences of 



pendences of the maximum positions for the square and 
triangular lattices are shown in the inset. At h/J < 1 
both C^(h/J) and C^(h/J) can be approximated by a 
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T/J 



FIG. 10: Temperature dependences of the heat capacity for 
HFM on a triangular lattice at h/J = 0.1, 0.5, 1 .0, 1.5 (from 
left to right). The present theory (solid), HTSEp 1 (dashed), 
and RPA (dotted). Thin lines correspond to h =0. 




12 3 

T/J 



FIG. 12: Temperature dependences of the heat capacity for 
HFM on a triangular lattice (from bottom to top): h/J =0.1, 
L =4, 6, 8, 10 (solid) and h/J =1, L =4, 6 (dashed). 



crease linearly 



Cm — 



0.51 



0.48 



0.45 - 



0.42 




FIG. 11: Field dependences of the heights and positions of 
the heat capacity maxima for HFM on square and triangular 
lattices. 



function 



with 



Cm — 



h 
J' 



0.5136, 
0.5254, 



0.0239, 
0.0189, 



square; 
triangular. 



A = 



-0.00693, 
-0.0248, 



B = 



B, 

0.5118, square; 
0.5468, triangular. 



Since decreases slower than C^, the inequality C^ > 
valid for low fields changes into the opposite one at 
h/J>2. 

Let us consider now the dependence of the thermody- 
namic functions on the cluster size L x L. It is inter- 
esting to determine the linear size Lq corresponding to 
the thermodynamic limit at a given magnetic field. This 
quantity is important, for example, on using such meth- 
ods as Monte Carlo and exact diagonalization, when a 
knowledge of an optimal cluster size makes it possible 
to obtain the thermodynamic functions of the infinite 
system within a reasonable volume of calculations. The 
dependence of the thermodynamic functions on L is also 
of practical interest, because of the isle structure of 3 He 
layers at some coverages.^ 

Fig. [12] displays temperature dependences of the heat 
capacity at h/J =0.1 and 1 for different cluster sizes L 
up to Lq. It is seen that with decrease in L maximum in 
the curve C(T) decreases and shifts to higher tempera- 
tures. At small L and very low fields a second maximum 
arises on the low temperature part of the heat capacity. 
A similar additional maximum resulting from the finite 
size of the system was found by ED for the 4x4 square 
lattice in Ref. [30j This result is also reproduced by our 
calculations. 

Fig. 



13 shows dependences L (h/J) for the square and 



At higher fields (h/J > 1.4) the maximum values de- triangular lattices. At h/J < 0.2 even small variation in 
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Our results agree well 




FIG. 13: Dependences Lo(h/J) for triangular (A) and square 
(■) lattices. 




FIG. 14: Temperature dependences of the magnetization for 
16 x 16 triangular lattice HFM at h/J =0.429, 0.214, 0.0429, 
0.0171, 6.09 10" 3 , 2.76 10" 3 , 4.29 10" 4 (from top to bottom): 
present theory (solid) and QMCp^l (symbols) . 



field leads to a sufficient change in Lq. As the field in- 
creases this dependence weakens, so that beginning with 
h/J ~ 0.2 rather small-sized clusters are appropriate for 
the numerical simulations of the real infinite systems. 
Temperature dependences of the magnetization for a 



16x16 triangular lattice together with the corresponding 
QMC datar^ are shown in Fig. 
with QMC. 

Now we check up the two criteria outlined in Sec. II, as 
applied to the triangular lattice HFM. The limiting value 
of the entropy is equal to 0.708 and 0.713 at h/J=0.05 
and 1, respectively, which slightly exceeds In 2. At low 
and high temperatures the function R(T, h) is close to 
zero as it was for HFM on the chain and square lat- 
tice. At fixed field in the intermediate temperature re- 
gion R(T, h) has a maximum, whose height decreases as 
h/J increases. The maximum value of R is ^0.094 at 
h/J = 0.05, whereas at h/J = 1 it does not exceed 0.046. 



IV. SUMMARY 

The thermodynamics of the low dimensional spin- 1/2 
Heisenberg ferromagnets in an external magnetic field is 
investigated within a second-order two-time Green func- 
tion formalism in the wide temperature and field range. 
The self-consistent set of equations for the correlation 
functions, vertex parameters and magnetization is ob- 
tained in the universal form appropriate for the descrip- 
tion of low dimensional HFM on a chain, square and 
triangular lattices. The fundamental point of our con- 
sideration is the account of the correct analytical prop- 
erties for the approximate transverse commutator Green 
function, from which the equation for the magnetization 
follows. This enables us to extend the range of adequate 
description for the HFM thermodynamics to lower fields 
as compared to the scheme proposed in Ref. [39] 

The thermodynamics of a triangular lattice HFM in 
a magnetic field is studied within a second-order Green 
function formalism for the first time. The temperature 
dependences of the magnetization, susceptibility, correla- 
tion functions, and heat capacity at different values of the 
magnetic field are calculated and analyzed in detail. For 
square and triangular lattices the positions and heights 
of the heat capacity maxima vs field are obtained. The 
dependences of the thermodynamic functions of the 2D 
HFM on the cluster size are investigated. For both types 
of lattices the cluster sizes corresponding to the thermo- 
dynamic limit are found as functions of field. 

The temperature and field dependences for the ther- 
modynamic functions calculated within our scheme are in 
close agreement with the corresponding results obtained 
by Bethe ansatz, quantum Monte Carlo simulations, high 
temperature series expansion, and exact diagonalization. 
Thus, we can conclude that the scheme used in this paper 
provides a good quantitative description for the thermo- 
dynamics of the low dimensional HFM in an external 
magnetic field on the three considered types of lattices 
for infinite as well as for finite-sized systems. 
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